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ABSTRACT 

We present constraints on the evolution of the virial to stellar mass ratio of galaxies 
with high stellar masses in the redshift range 0.2 < z < 0.8, by comparing weak lensing 
measurements of virial mass Mvir to estimates of stellar mass Afstar from COMBO-17. For a 
complete sample of galaxies with log(A/stai /A^0) > 10.5, where the majority show an early- 
type morphology, we find that the virial mass to stellar mass ratio is given by Mvir/Mgtar — 
53^\q. Assuming a baryon fraction from the concordance cosmology, this corresponds to 
a stellar fraction of baryons in massive galaxies of 57^/57b = 0.10 ± 0.03. Analysing the 
galaxy sample in different redshift slices, we find little or no evolution in the virial to stellar 
mass ratio, and place an upper limit of ~ 2.5 on the growth of massive galaxies through the 
conversion of gas into stars from z = 0.8 to the present day. 



1 INTRODUCTION 

In the standard cold dark matter paradigm, luminous galaxies re- 
side in massive dark matter halos. The large-scale clustering and 
formation of these luminous galaxies is driven by the clustering 
and formation histories of their dark matter halos. In contrast, the 
visible properties of such galaxies are a strong function of the com- 
plex physics of baryonic matter, in particular the laws governing 
stellar and AGN feedback, gas cooling and star formation. The lu- 
minosities and colors of galaxies are a strong function of the ef- 
ficiencies of these processes, and our inability to robustly model 
them is a critical limitation of current models of galaxy formation 
and evolution (as discussed, for example, bv Benson et aLi.2003il . 
A key way to constrain the efficiency of such processes is to mea- 
sure empirically the ratio of stellar mass to the total virial mass 
of the halo and galaxy system, for certain subsets of the galaxy 
population. Assuming a baryon to dark matter fraction from the 
concordant cosmology, it is then possible to determine empirically 
the fraction of baryons which manage to cool and condense into 
stars in halos of a given mass. Such measurements of the mass out 



to ^ 100 kpc scales not only constrain the efficiency of complex 
'gastrophysical' processes, but they also empirically calibrate and 
test halo oc cupation distribut ion analyses of galaxy properties (see 
for example lYang et all200^. 

Satellite galaxies jzaritskv & Whitd 1 19941: IPrada et alJ 

[200j) and weak gravitation al lensing (see the review by 
iBartelmarm & Sch neider 2001"), are currently the most promising 
methods available to probe virial masses of dark matter halos. 
In this letter we focus on weak gravitational lensing, the weak 
distortion of distant sources by the gravitational field of foreground 
structure, as it provides a unique way to probe matter on all 
scales, at all redshifts, irrespective of its dynamical state or nature. 
Termed 'galaxy-galaxy lensing', the weak tangential shearing of 
background galaxies relative to foreground galaxies can therefore 
be used to probe the distribution of dark matter in galaxies, as first 
demonstrated by Brainerd et al. ( 1996). 

In this letter we combine information from two surveys of 
the Ch andra Deep Field South; the ground-based C OMBO-17 
surve y iWolf et aljEool and the HST GEMS survey <Rix et alJ 
|2004) . in order to constrain the virial to stellar mass ratio in stel- 
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lar massive galaxies out to z; ~ 0.8. The 17-filter photometry of 
COMBO- 17 provides accurate photometric redshifts with errors 
cjz ~ 0.02(1 + z). COMBO-17 additionally provides stellar mass 
estimates Matar by fitting a parameterised star formation history 
to the lo w resolution 17-band galaxy spectra iBorch et al. 2006), 
based on a iKroupa etalJ<1993h initial mass function (Kroupa IMF). 
The HST resolution of the GEMS imaging allows for high signal- 
to-noise measurements of the 'galaxy-galaxy' lensing signal which 
can provide an ensemble estimate of the virial mass for (sub-)sets of 
galaxies. A detailed description of the galaxy-galaxy lensing anal- 
ysis of GEMS and the resulting constraints on different dark matter 
halo profiles for different galaxy types will be presented in a forth- 
coming paper (Heymans et al. in prep). For the purposes of this 
letter however we focus on a complete sample of 626 lens galaxies 
with high stellar masses. We measure the virial to stellar mass ra- 
tio as a function of galaxy redshift to z ~ 0.8, presenting the first 
redshift dependent galaxy-galaxy lensing analysis to be performed 
consistently over such a large redshift range. 

This letter is organised as follows. In Section|2|we describe the 
theory that underpins our measurement of virial masses and review 
the particulars of our analysis. We refer the reader to section 3 of 
[Heymans et al. 1 2005) hereafter H05, for a lensing specific descrip- 
tion of the GEMS observations and the data reduction. In Section|3| 
we present constraints on the virial radius and the virial to stellar 
mass ratio for our full stellar massive galaxy sample and convert 
this into a stellar fraction of baryons. We compare our results to the 
low-redshift re sults of .Hoeks tra et al. 1 2005), hereafter HHYLG, 
andlMandelba um et aljl2006l) . hereafter M06, and investigate the 
dependence of the virial to stellar mass ratio on galaxy redshift, 
setting constraints on the growth of massive galaxies through the 
conversion of gas into stars from z ~ 0.8 to the present day. 
We discuss the implications of these results and conclude in sec- 
tion|4] Throughout this letter we assume a ACDM cosmology with 
Q.rn. = 0.3, f^A = 0.7 and Ho = 100/ikms~^ Mpc~\ For our 
virial to stellar mass ratio measurement we set = 0.7. 



2 MEASURING VIRIAL MASSES WITH 
GALAXY-GALAXY LENSING 

The mass of a foreground Tens' galaxy will weakly shear the im- 
ages of background or 'source' galaxies. The extent of this tan- 
gential shearing as a function of lens and source position can be 
directly derived from the density distribution of the lens galaxy. In 
this letter we adopt the NEW profile as a model for dark matter ha- 
les, where the density distribution of a halo at redshift z is given 
by 

pM = ^-^P^ (1) 

5c is the characteristic density, is the scale radius and pc{z) is 
the critical density^ given by 3H(z)^ /8tvG. The virial radius is 
defined as the radi us where the mass density of the halo is equal to 
A X o^. lEke et alj n996) show that for flat A cosmology, assuming 
a spherical collapse model, 

A{nm,A) = 178nV'^ . (2) 



Note that within the galaxy-galaxy lensing literature pc is often defined 
differently. In this le tter we follow the definition given in the Appendix of 
iNavarro et ai]<199ii . 



For the cosmology adopted here A « 100 at 2 = 0, and A ~ 150 
at z — 0.8. The corresponding virial mass'^ Ma is given by 

Ma = Apc^rl . (3) 

Defining the concentration parameter as c = r/s/rs, the character- 
istic halo density Sc is then given by 

^ A 

" ^ 3 ln(l + c) - c/(l + c) • ^ ' 

Note that for a given CDM cosmology, halo mass Myir and con- 
centration c, or scale radius r^, ar e related jNavarroeta'i]ll997t 
iBullock et all200lLlEke et alj200ll) . where the dependence is cal- 
culated through fits to numerical simulations. 

The expression for the wea k lensing shear 7 induced by an 
NFW dark matter halo is given in IWright & Brainerjj200(]l) . 

2.1 Lens galaxy sample selection 

To measure accurate virial masses it is important to have accurate 
photometric redshifts. To compare virial mass to stellar mass we 
also require good stellar mass estimates. We therefore limit our- 
selves to a sample of lens galaxies with a measured COMBO-17 
photometric redshift (this implies a magnitude limit of i? < 24) 
between 0.2 < z < 0.8. For the purposes of this letter we wished 
to create a complete sample of galaxies in order to be able to study 
the virial to stellar mass ratio as a function of redshift without 
the inclusion of potentially biasing selection effects. This require- 
ment imposed a high stellar mass cut where log(Mstar/M0) > 
10.5, above which it has been shown that the COMB O-17 cata- 
logue is complete in both redshift estimates and flux ^W olf etalJ 
[2004; Borch et al. 2006). The Sersic fits of these galaxies (see 
Barden et al. 2005; Mcintosh et al. 2005) indicate that the majority 
of the lens sample (~ 75%) are early-type galaxies with n > 2.5. 
Note that including an n > 2.5 Sersic selection criteria such that 
our lens sample contains only early-type galaxies does not change 
our results except for decreasing the lensing signal-to-noise. The 
galaxy sample has an average photometric redshift accuracy of 
~ 0.03 and a mean stellar mass of (Afatar) = 7.2 X W^°Mq. 
The galaxies are very luminous with a mean SDSS r-band lumi- 
nosity (L) = 2.4 i*, where i* = lO^°/i"^L0. They also have a 
low dynamic range of galaxy luminosity with o-(log L/L*) — 0.2. 

Source galaxies with which to measure the galaxy-galaxy 
lensing signal are selected from the GEMS data as described in sec- 
tion 4.3 of H05. This yields a number density of 65 source galaxies 
per square arcmin. 

The lens redshift limits were chosen to optimise the measured 
signal. Galaxies at higher redshifts typically have poorer photomet- 
ric redshift accuracy and they also provide rather weak constraints 
owing to the reduced number of background sources. The lower 
redshift limit results from the fact that we exclude all sources for 
which lenses might lie outwith the GEMS area. For a given max- 
imum impact parameter rmax, this corresponds to a minimum an- 
gular separation 0{zi^in) from the edge of the field. We therefore 
include the lower redshift limit to reduce 9{zi min) thus including 
as many sources as possible. 

In this letter we will make two assumptions. As our lens sam- 
ple spans a low dynamic range in luminosities we will assume that 

^ This definition of virial mass is directly equivalent to that used by 
HHYLG as A{z)pc{z) = Avir(z)pbg(z) where Avir and p^g are given 
after equation 9 of HHYLG. 
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the halo virial radius r a of our lens sample is constant. M06 find 
that stellar mass is a good galaxy-type-dependent proxy for halo 
mass. We will therefore also assume that, on average, the halo virial 
mass of our galaxy sample scales linearly with galaxy stellar mass 
by defining the following scaling relation, 



Mvir = fc^ 



/Ms, 



(5) 



Using the mass/concentration relations from lNavarro et all il997h . 
the halo profile of a galaxy in our sample can then be uniquely 
described by an NFW model with one free parameter; either r a or 

fcNPW. 



2.2 Analysis Method 

We u se the metho d of iKaiser et alj (^^|), iLuDpino & KaiseJ 
(l9^ and Hoekstr a~alJ Jl998h (KSB-t-l to obtain an estimate of 

the observed ellipticity 6°'^'' of each source galaxy in the GEMS sur- 
vey. In the weak shear limit, where 7^1, the observed ellipticity 
6°*^^ ^ e° +7, where is the intrinsic source galaxy ellipticity. The 
application of the method to the GEMS data set is described in de- 
tail in H05 along with the results of several diagnostic tests which 
confirm the success of this method at removing PSF distortions. 
The levels of shear calibration bias expected with t he HQS imple- 
menta tion of KSB-l- has been shown to be ~ 3% I Hevmans et alJ 
l2006h which is well within the statistical uncertainties of this anal- 
ysis. 

In order to interpret the measured average tangential shear 
in the context of NFW halo lensing, redshift information for 
bot h the lens and source g alaxies is required (see for exam- 
ple |Hei^ei^c^ni]|20Ql). The COMBO- 17 survey provides 
accurate photometric redshift s for the foreg round lens sample 
of galaxies with i? < 24 JWoIf et alj '2004'). For the fainter 
source galaxy sample only the redshift probability distribution is 
known (see section 6 of H05). The maximum-likelihood method of 
ISchneider & Rixl il997h was designed to take advantage of such a 
data set by analysing the galaxy-galaxy lensed shear statistically, 
and it is this method that we use in the following analysis and de- 
scribe below. 

For a model galaxy density profile, in the case where all galaxy 
redshifts are known, the weak shear 7 experienced by each source 
galaxy can be predicted by summing up the shear contributions 
from all the foreground lens galaxies, taking into account the multi- 
ple weak deflections that a source may experience liBrainerd.2005.V 
Note that when a sample of lenses are selected from the foreground 
galaxy population, the other foreground galaxies act as a source 
of noise without introducing bias ( Kleinheinrich et al. 2006). In 
this analysis the redshifts of the source galaxies are unknown, and 
we therefore assign those galaxies a magnitude-dependent redshift 
probability distribution p(2,mag) (equation 15 H05) and calcu- 
late the expectation value of the shear (7). This is done through 
Monte Carlo integration, drawing a source galaxy redshift estimate 
z'^ from the distribution p[z,m&g), v = L.A^mc times, where 
A'mc = 100 in this analysis. For each estimate the induced 
galaxy-galaxy lensing shear 7" is calculated with the resulting 
mean shear given by 



(7> = 



iVi 



MC 



(6) 



Note that in practice we only consider the shear from those lenses 
where the angular separation from the source Q < rmax/ Di with 
fmax = 300ft~^kpc. The intrinsic source galaxy ellipticity is 
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Figure 1. The average tangential alignment (74) of source galaxies with 
respect to the lens sample as a function of angular separation on the sky 
9. The solid curve shows the theoretical lensing signal from an NFW halo 
at z = 0.65, the median redshift of the lens sample, with the maximum 
likelihood method constrained virial radius ta = 204/i~^kpc. 



then calculated, £^ ~ e — 7. The distribution of each component 
of the observed galaxy ellipticity is well described, for the GEMS 
survey, by a Gaussian of width o-j = 0.31. As the induced shear 7 
is weak, the probability for observing an intrinsic ellipticity of is 
then given by 



27r cr; 



■ exp 



2cr| 



(V) 



The best-fit dark matter halo parameters are determined by max- 
imising the likelihood L = 11 [P(e^)i] where the product extends 
over all source galaxies i. 



3 RESULTS 

In this section we present our constraints on the virial radius and 
the virial to stellar mass ratio of the stellar massive galaxies in our 
lens sample. Assuming that the chosen subset of lens galaxies have 
a constant virial radius we find ta = 2Q4^22'''^^kpc. This cor- 
responds to a virial mass of Mvir = 29.8lg g lO^^/i~^M0 for a 
galaxy at z ~ 0.65, the median redshift of the lens sample. Note 
that the corresponding virial mass is dependent on the redshift of 
the halo because of the redshift dependence of the critical density 
pc and A (equation |3}- For our mean galaxy luminosity we find 
an average mass-to-light ratio for the stellar massive galaxies of 
Mvir/i = 123 ± 36/1M0/L0. To illustrate this result, obtained 
using the maximum likelihood method described in section lT2l fig- 
ureQshows the average tangential alignment {-yt) of source galax- 
ies with respect to the lens sample as a function of angular separa- 
tion on the sky. The measured signal (74) is well fit by the expected 
signal around an NFW halo at the median redshift of the lens sam- 
ple with the maximum likelihood method constrained virial radius 
of r-A = 204/i"^kpc. 

Our results are in good agreement with the results of HHYLG 
who find Mvir = 23.3l5;i x lO"/i~^M0 for their luminous 
galaxy sample at z ~ 0.32, where (L) = 2.491/,. Scaling our 
results to 2 = 0.32 we find Mvir = '^T-^tfJj x lO^/i^^Mg. 
These results are also in good agreement with the results of M06 
who find a virial mass for central galaxies (as opposed to satellite 



galaxies) of Ma 



14.ll^ j lO^^/i-^M0, for a sample of early- 



type galaxies with mean stellar mass, (Msi 



5.8 X lO^^Mp 
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and mean luminosity (L) ~ l.SL*. This galaxy sample has a 
mean redshift z ~ 0.11, and scaling our results to this redshift 
we find a consistent virial mass of Mvir = 12.5t4 2lO"/i"^Mo 
for our sample of stellar massive galaxies with mean stellar mass, 
(M,tar) = 7.2 X 1O'°M0 and mean luminosity (L) ~ 2.4L*. 

In Section im we introduced a stellar mass scaling to relate 
halo mass to stellar mass, which we choose to be linear"^. For the 
full lens sample we measure a best fit value for the NFW virial 
mass to stellar mass ratio as fcNPw ~ SSlJ^g with Icr errors. Note 
that AiNFW = 1 is ruled out with 99.99% confidence. Our result is 
in good agreement with the results of HHYLG and M06, as shown 
in Figure |2| HHYLG find fcNPW ~ 45 ± 6 for a scaled Salpeter 
IMF which is similar to the Kroupa IMF used in this analysis'*. 
M06 find fcNFW = 35 ± 13 using a Kroupa IMF, for the sample 
of early-type galaxies whose average stellar mass best matches the 
average stellar mass of our sample, as discussed above. For a higher 
stellar mass selected sample with (A/star) = 11-2 x lO"'!/©, M06 
find A:nfw = 43 ± 13. As ~ 20% of our sample would fall into 
this mass bin, we also show this result in Figure|2| 

Following the analysis of HHYLG and M06 we convert our 
virial to stellar mass ratio into a stellar baryon fraction /, by as- 
suming that the baryon to dark matter ratio in massi ve galaxies is 
given by the global value ^h/^m ~ 0.176 ± 0.013 iSoergel et alj 
l2006t) such that 



m 



0,1/0,1, 



(8) 



Mvir Ob 

The measured virial to stellar mass ratio then gives us a stel- 
lar baryon fraction of /« = 01 /Ob — 0.10 ± 0.03, suggest- 
ing that galaxies with high stellar masses are rather inefficient at 
converting baryons into stellar mass: only 10% of the baryons 
are contained in stars. T his result is in good a greement with the 
cosmic baryon budget of Fuk ugita et al r il99^ ^ who estimate the 
the amount of baryons in different states, finding the fraction of 
baryons in spheroids to be /» ~ 0.12. 



3.1 Redshift dependence 

Splitting our galaxy sample into four redshift slices we show, in 
Figure|2| that the virial to stellar mass ratio remains fairly constant 
as a function of galaxy redshift. For a linear fit to these results, 
together with the estimates of HHYLG and M06, we can constrain 
kNFwiz) = (34± 11) + (31±35)z. Comparing the virial to stellar 
mass ratio at z = and z — 0.8, we find a la confidence region 
where 0.8 < {k[z = 0.8]/k[z = 0]) < 2.6. The data therefore 
indicate that the virial to stellar mass ratio has evolved by a factor 
< 2.5 since z=0.8. 



4 DISCUSSION AND CONCLUSION 

In this letter, we use a weak gravitational lensing study of the HST 
GEMS survey, in combination with the 17-colour photometry of the 
COMBO- 17 survey, to determine the ratio between virial and stel- 
lar mass in a complete sample of massive galaxies out to 2: ~ 0.8. 

^ Note that we also tested non-linear stellar mass scaling, but found for 
our complete sample of log(Mstar/Af0) > 10.5 galaxies, any non-linear 
component was consistent with zero within the noise 
^ Note that the Kroupa IMF and scaled Salpeter IMF are roughly half the 
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Figure 2. The virial mass to stellar mass ratio as a function of galaxy red- 
shift (circles). The result at z = 0.32 is from HHYLG (star). The results at 
z = 0.11 and z = 0.13 are from M06 (triangle) measured from different 
stellar mass selected samples that best match our lens sample. The errors are 
1(T in all cases, with the exception M06 who quote 3cr errors. The hatched 
region is the Icr confidence regions for the virial mass to stellar mass ratio 
measured from GEMS over the full redshift range 0.2 < z < 0.8. The 
solid line indicates the best linear fit to Mvir/Mstar(z) and the dashed 
lines denote the Icr eiTors to this fit. 



We find a mean virial to stellar mass ratio A^vir/A^star — 53^J5 
with little or no redshift evolution in this stellar mass fraction. 

This result provide s a high redshi ft extension to the analysis 
of iHoekstra etai] J2005l) (HHYLG) and lMandelbaum et alj Ew^ 
(M06). The excellent agreement found between the results of this 
analysis, HHYLG and M06 is important to comment upon as the 
analyses were performed using very different techniques and came 
from very different surveys: HHYLG use the Red sequence Clus- 
ter Survey, a medium-deep 45.5 square degree ground-based sur- 
vey in BVRz, M06 use the shallow Sloan Digital all Sky Sur- 
vey which is ground-based in ugriz. This analysis uses GEMS, 
a deep 0.25 square degree space-based survey with accompanying 
17-band COMBO- 17 data. 

In this analysis we use stellar masses calculated from the best- 
fitting models of star formation history to the low resolution 17- 



comparison to the HHYLG Salpeter IMF measurement of fci; 



: 26±4. 



band galaxy spectra. M06 use a similar technique 1 Kauffmarir^^I 
|20^). HHYLG use the simple relationship of iBell & de Jong 
(2001) that relates stellar mass to B ~V colour. This relationship 
has a systematic uncertainty that is often estimated to be ~ 30% 
(.Bell et aL 2003 ). Stellar mass uncertainty is indeed hard to avoid 
with all methods, however, owing to the unknown f requency of 
bursts of star formation and uncertainty in dust models iBorch et alJ 
I2OO6). The fact all three results are in such close agreement, how- 
ever, demonstrates that the simple models of HHYLG are con- 
sistent with the more complex star formation history models of 
COMBO-17 and M06, which bodes well for similar studies in fu- 
ture broad-band surveys. 

We have applied the lSchneider & Rixlll997i) maximum likeli- 
hood method to set constraints on the virial radius and the virial to 
stellar mass ratio using an NFW dark matter profile. This method 
takes into account the uncertainty in the redshift of both the source 
and lens galaxies, using a Monte Carlo technique and includes the 
effects of the multiple deflection of light by successive foreground 
lens galaxies. In contrast, HHYLG measure the tangential shear 
around foreground galaxies for different luminosity bins and fit this 
result with the shear expected from an NFW profile. The HHYLG 
galaxies are selected to be isolated such that the lensed background 
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galaxies experience only the gravitational potential of the isolated 
galaxy. Their larger photometric redshift errors are more likely to 
scatter faint galaxies into higher luminosity bins than scatter bright 
galaxies into lower luminosity bins. This then leads to a bias in 
the measurement of the mass at fixed luminosity to a lower value. 
HHYLG therefore need to apply a correction factor to the virial 
masses which they determine from an analysis of mock catalogues. 
This bias is not a problem for our analysis owing to the higher ac- 
curacy of photometric redshifts and the use of the Schneide r & Rixl 
jl997|) method. Our observational systematic errors are therefore 
minimal, and we are limited by random errors from smaller num- 
ber statistics. The third analysis method of M06 uses the halo 
model to extract information from the galaxy-galaxy lensing sig- 
nal and separate contribu tions from central and satellite galaxies 
jMandelbaum et alj|2005l) . The fact that these three very different 
methods produce consistent results serves as a good indication that 
the three results are robust. 

This analysis is the first time a comparison of the properties of 
luminous galaxies and their parent dark matter halos has been per- 
formed consistently across a substantial span of cosmic time. The 
galaxy sample has been selected to be complete out to z ~ 0.8 
and we are therefore not subject to redshift-dependent selection 
biases. We have shown that for this galaxy sample, the majority 
of which show early-type morphologies, the virial to stellar mass 
ratio as a function of galaxy redshift remains fairly constant with 
M„ir/Mstar ~ 50: Specifically, we place a Icr upper limit on the 
growth of the stellar to virial mass ratio since z = 0.8, of a factor 
< 2.6. Because we probe the ratio of stellar to virial mass, we can- 
not probe merger-driven growth of galaxies as stellar masses and 
virial masses grow in lock-step. Instead, our constraint effectively 
limits the conversion of gas into stars in the massive galaxy pop- 
ulation, i.e. star formation, as a function of redshift. We therefore 
can conclude that the growth of stellar mass due to star formation 
is limited, to a factor of < 2.5 from z = 0.8 to the present day. 
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